Genome Assembly Research Articles

Understanding bacterial pathogen diversity: A proteogenomic analysis and use of an array of genome assemblies to identify novel virulence factors of the honey bee bacterial pathogen Paenibacillus larvae.

Mass spectrometry proteomics data are typically evaluated against publicly available annotated sequences, but the proteogenomics approach is a useful alternative. A single genome is commonly utilized in custom proteomic and proteogenomic data analysis. We pose the question of whether utilizing numerous different genome assemblies in a search database would be beneficial. We reanalyzed raw data from the exoprotein fraction of four reference Enterobacterial Repetitive Intergenic Consensus (ERIC) I-IV genotypes of the honey bee bacterial pathogen Paenibacillus larvae and evaluated them against three reference databases (from NCBI-protein, RefSeq, and UniProt) together with an array of protein sequences generated by six-frame direct translation of 15 genome assemblies from GenBank. The wide search yielded 453 protein hits/groups, which UpSet analysis categorized into 50 groups based on the success of protein identification by the 18 database components. Nine hits that were not identified by a unique peptide were not considered for marker selection, which discarded the only protein that was not identified by the reference databases. We propose that the variability in successful identifications between genome assemblies is useful for marker mining. The results suggest that various strains of P. larvae can exhibit specific traits that set them apart from the established genotypes ERIC I-V.

Distribution and expression of the aac(6')-Im (aacA16) aminoglycoside resistance gene.

The aac(6')-Im (aacA16) amikacin, netilmicin and tobramycin resistance gene cassette had been circulating globally undetected for many years in a sublineage of Acinetobacter baumannii global clone 2. To identify sources for the aac(6')-Im fragment found in A. baumannii. MinION long-read sequencing and Unicycler hybrid assemblies were used to determine the genetic context of the aac(6')-Im gene. Quantitative reverse transcriptase PCR was used to measure expression. Among >60 000 non-Acinetobacter draft genomes in the MRSN collection, the aac(6')-Im gene was detected in Pseudomonas putida and Enterobacter hormaechei isolates recovered from patients in Thailand between 2016 and 2019. Genomes of multiply resistant P. putida MRSN365855 and E. hormaechei MRSN791417 were completed. The class 1 integron containing the aac(6')-Im cassette was in the chromosome in MRSN365855, and in an HI2 plasmid in MRSN791417. However, MRSN791417 was amikacin susceptible and the gene was not expressed due to loss of the Pc promoter of the integron. Further examples of aac(6')-Im in plasmids from or the chromosome of various Gram-negative species were found in the GenBank nucleotide database. The aac(6')-Im context in integrons in pMRSN791417-8 and a Klebsiella plasmid pAMR200031 shared similarities with the aac(6')-Im region of AbGRI2-Im islands in A. baumannii. In other cases, the cassette array including the aac(6')-Im cassette was different. The aac(6')-Im gene is widespread, being found so far in several different species and in several different gene cassette arrays. The lack of amikacin resistance in E. hormaechei highlights the importance of correlating resistance gene content and antibiotic resistance phenotype.

Complete circular assembly of Histophilus somni PDS1537697-81-1 genome and phenotypic characterization of its antibiotic susceptibility.

A Histophilus somni isolate from a clinically healthy, fall-placed calf was obtained upon arrival to a commercial feedlot. Fall-placed calves are commonly viewed to be at high risk for the development of bovine respiratory disease. The isolate was phenotyped for antimicrobial susceptibility and sequenced to obtain a complete, circular, genome assembly.

Different control of resistance to two Colletotrichum orbiculare pathogenic races 0 and 1 in cucumber (Cucumis sativus L.).

Quantitative trait locus analysis identified independent novel loci in cucumbers responsible for resistance to races 0 and 1 of the anthracnose fungal pathogen Colletotrichum orbiculare. Cucumbers have been reported to be vulnerable to Colletotrichum orbiculare, causing anthracnose disease with significant yield loss under favorable conditions. The deployment of a single recessive Cssgr gene in cucumber breeding for anthracnose resistance was effective until a recent report on high-virulent strains infecting cucumbers in Japan conquering the resistance. QTL mapping was conducted to identify the resistance loci in the cucumber accession Ban Kyuri (G100) against C. orbiculare strains 104-T and CcM-1 of pathogenic races 0 and 1, respectively. A single dominant locus An5 was detected in the disease resistance hotspot on chromosome 5 for resistance to 104-T. Resistance to CcM-1 was governed by three loci with additive effects located on chromosomes 2 (An2) and 1 (An1.1 and An1.2). Molecular markers were developed based on variant calling between the corresponding QTL regions in the de novo assembly of the G100 genome and the publicly available cucumber genomes. Multiple backcrossed populations were deployed to fine-map An5 locus and narrow the region to approximately 222 kbp. Accumulation of An2 and An1.1 alleles displayed an adequate resistance to CcM-1 strain. This study provides functional molecular markers for pyramiding resistance loci that confer sufficient resistance against anthracnose in cucumbers.

Genomic Insights Into Enterococcus mundtii 203: A Promising Probiotic Candidate Isolated From Camel Feces.

Enterococcus, a common commensal organism in the human gut, exhibits a dual nature with certain strains offering probiotic benefits, while others are associated with nosocomial infections. In this study, we conducted a comprehensive examination of the genome of Enterococcus mundtii strain 203 to assess its probiotic potential and safety profile. The complete genome sequencing, assembly, and annotation were performed, followed by bioinformatics analysis. Our investigation reveals a detailed characterization of the Enterococcus mundtii 203 genome, originally isolated from camel feces, with a size of 3,053,234 bases and a GC content of 38.4%. Importantly, our analysis suggests that this strain poses no risk as a human pathogen due to the absence of antibiotic resistance determinants and virulence factors. The genome harbors a multitude of genes responsible for lactic acid production, bioactive peptide synthesis, adhesion molecule expression, resistance to harsh gut conditions, and enhancement of host metabolism. These findings underline the potential probiotic functionality of Enterococcus mundtii 203, positioning it as a promising candidate. Notably, our study did not identify any sequences related to insertion elements or CRISPR-Cas fragments.

Comparative genomic analyses provide insight into the pathogenicity of three Pseudomonas syringae pv. actinidiae strains from Anhui Province, China

BackgroundPseudomonas syringae pv. actinidiae (Psa) is an important bacterial plant pathogen that causes severe damage to the kiwifruit industry worldwide. Three Psa strains were recently obtained from different kiwifruit orchards in Anhui Province, China. The present study mainly focused on the variations in virulence and genome characteristics of these strains based on the pathogenicity assays and comparative genomic analyses.ResultsThree strains were identified as biovar 3 (Psa3), along with strain QSY6 showing higher virulence than JZY2 and YXH1 in pathogenicity assays. The whole genome assembly revealed that each of the three strains had a circular chromosome and a complete plasmid. The chromosome sizes ranged from 6.5 to 6.6 Mb with a GC content of approximately 58.39 to 58.46%, and a predicted number of protein-coding sequences ranging from 5,884 to 6,019. The three strains clustered tightly with 8 Psa3 reference strains in terms of average nucleotide identity (ANI), whole-genome-based phylogenetic analysis, and pangenome analysis, while they were evolutionarily distinct from other biovars (Psa1 and Psa5). Variations were observed in the repertoire of effectors of the type III secretion system among all 15 strains. Moreover, synteny analysis of the three sequenced strains revealed eight genomic regions containing 308 genes exclusively present in the highly virulent strain QSY6. Further investigation of these genes showed that 16 virulence-related genes highlight several key factors, such as effector delivery systems (type III secretion systems) and adherence (type IV pilus), which might be crucial for the virulence of QSY6.ConclusionThree Psa strains were identified and showed variant virulence in kiwifruit plant. Complete genome sequences and comparative genomic analyses further provided a theoretical basis for the potential pathogenic factors responsible for kiwifruit bacterial canker.

Detection of Anaplasma and Ehrlichia bacteria in humans, wildlife, and ticks in the Amazon rainforest

Tick-borne bacteria of the genera Ehrlichia and Anaplasma cause several emerging human infectious diseases worldwide. In this study, we conduct an extensive survey for Ehrlichia and Anaplasma infections in the rainforests of the Amazon biome of French Guiana. Through molecular genetics and metagenomics reconstruction, we observe a high indigenous biodiversity of infections circulating among humans, wildlife, and ticks inhabiting these ecosystems. Molecular typing identifies these infections as highly endemic, with a majority of new strains and putative species specific to French Guiana. They are detected in unusual rainforest wild animals, suggesting they have distinctive sylvatic transmission cycles. They also present potential health hazards, as revealed by the detection of Candidatus Anaplasma sparouinense in human red blood cells and that of a new close relative of the human pathogen Ehrlichia ewingii, Candidatus Ehrlichia cajennense, in the tick species that most frequently bite humans in South America. The genome assembly of three new putative species obtained from human, sloth, and tick metagenomes further reveals the presence of major homologs of Ehrlichia and Anaplasma virulence factors. These observations converge to classify health hazards associated with Ehrlichia and Anaplasma infections in the Amazon biome as distinct from those in the Northern Hemisphere.

An integrated multi-omics approach reveals polymethoxylated flavonoid biosynthesis in Citrus reticulata cv. Chachiensis

Citrus reticulata cv. Chachiensis (CRC) is an important medicinal plant, its dried mature peels named “Guangchenpi”, has been used as a traditional Chinese medicine to treat cough, indigestion, and lung diseases for several hundred years. However, the biosynthesis of the crucial natural products polymethoxylated flavonoids (PMFs) in CRC remains unclear. Here, we report a chromosome-scale genome assembly of CRC with the size of 314.96 Mb and a contig N50 of 16.22 Mb. Using multi-omics resources, we discover a putative caffeic acid O-methyltransferase (CcOMT1) that can transfer a methyl group to the 3-hydroxyl of natsudaidain to form 3,5,6,7,8,3’,4’-heptamethoxyflavone (HPMF). Based on transient overexpression and virus-induced gene silencing experiments, we propose that CcOMT1 is a candidate enzyme in HPMF biosynthesis. In addition, a potential gene regulatory network associated with PMF biosynthesis is identified. This study provides insights into PMF biosynthesis and may assist future research on mining genes for the biosynthesis of plant-based medicines.

Aroma formation in single- and double-petal jasmines (Jasminum sambac) during flowering via volatile metabolome and transcriptome

Jasminum sambac (L.) Aiton exists in single-petal (SP), double-petal (DP), and multi-petal (MP) phenotypes, which have a sweet floral scent and release aroma when the flowers open. A total of 103 volatiles were discovered from SP and DP jasmines. 16 volatiles were identified as their key odorants. Caryophyllene, geraniol, linalool, α-farnesene, β-ocimene, methyl anthranilate, methyl salicylate, benzeneacetaldehyde and hexanal may form the overall aroma of those two jasmines during flowering. High-resolution transcriptional profiles of those two jasmines at four flowering stages were generated from our recently reported chromosome-level genome assembly, from which we identified 35 differentially expressed genes (DEGs) involved in terpenoid metabolic pathways and 10 DEGs involved in the phenylpropane/benzene metabolic pathway, which may shape the characteristic aroma of those two jasmines at the molecular level. Through co-expression network analysis of key odorants, we found that bZIP8, bHLH6, MYB4, and AP2/ERF3 may contribute to the regulation of aroma volatiles in those two jasmines. These results offer an important theoretical foundation for studying the regulation of aroma formation for SP and DP jasmines.

The first high-quality genome assembly and annotation of Lantana camara, an important ornamental plant and a major invasive species

This study presents the first annotated, haplotype-resolved, chromosome-scale genome of Lantana camara, a flowering shrub native to Central America and known for its dual role as an ornamental plant and an invasive species. Despite its widespread cultivation and ecological impact, the lack of a high-quality genome has hindered the investigation of traits of both ornamental and invasive. This research bridges the gap in genomic resources for L. camara, which is crucial for both ornamental breeding programs and invasive species management. Whole-genome and transcriptome sequencing were utilized to elucidate the genetic complexity of a diploid L. camara breeding line UF-T48. The genome was assembled de novo using HiFi and Hi-C reads, resulting in two phased genome assemblies with high Benchmarking Universal Single-Copy Orthologs (BUSCO) scores of 97.7%, indicating their quality. All 22 chromosomes were assembled with pseudochromosomes averaging 117 Mb. The assemblies revealed 29 telomeres and an extensive presence of repetitive sequences, primarily long terminal repeat transposable elements. The genome annotation identified 83,775 protein-coding genes, with 83% functionally annotated. In particular, the study mapped 42 anthocyanin and carotenoid candidate gene clusters and 12 herbicide target genes to the assembly, identifying 38 genes spread across the genome that are integral to flower color development and 53 genes for herbicide targeting in L. camara. This comprehensive genomic study not only enhances the understanding of L. camara’s genetic makeup but also sets a precedent for genomic research in the Verbenaceae family, offering a foundation for future studies in plant genetics, conservation, and breeding.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator)

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as “avivorous bats”). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph.

Despite advances in long-read sequencing technologies, constructing a near telomere-to-telomere assembly is still computationally demanding. Here we present hifiasm (UL), an efficient de novo assembly algorithm combining multiple sequencing technologies to scale up population-wide near telomere-to-telomere assemblies. Applied to 22 human and two plant genomes, our algorithm produces better diploid assemblies at a cost of an order of magnitude lower than existing methods, and it also works with polyploid genomes.

An assembly sequence monitoring method based on workflow modeling for human–robot collaborative assembly

An assembly sequence monitoring method based on workflow modeling for human–robot collaborative assembly

Genomic Advancement in Wheat (Triticum aestivum L.): Harnessing Technological Breakthroughs for Future Strategies

Wheat can be greatly improved by identifying genes that are significant to agronomy. Although progress in wheat genetics and genomics has been limited by the genome’s vast size and complexity. Wheat has a high-quality genome sequence in light of recent developments in genome sequencing and sequence assembly. Here, we propose that wheat biology can benefit from the same approaches that have been used to define biological systems in model species, such as the generation and characterization of mutants, methods for gene cloning, and enhanced transgenic technology. These tactics will encourage the establishment of wheat breeding programs and hasten the field of wheat biology. We also summarize current developments in functional genomics of wheat. In order to contribute to global food security, we conclude by talking about the future of wheat functional genomics and the sensible design-based molecular breeding of new wheat varieties. We suggest that researchers studying wheat should embrace the methods utilized for functional genomic analysis in other model species. One method for doing this is using gene cloning to find the gene causing an intriguing phenotype, deciphering the biological role of genes through the analysis of their corresponding mutants, and creating a related wheat mutant from the proposed wild-type gene to confirm the target gene’s functionality.

Whole Snake Genomes from Eighteen Families of Snakes (Serpentes: Caenophidia) and Their Applications to Systematics.

We present genome assemblies for 18 snake species representing 18 families (Serpentes: Caenophidia): Acrochordus granulatus, Aparallactus werneri, Boaedon fuliginosus, Calamaria suluensis, Cerberus rynchops, Grayia smithii, Imantodes cenchoa, Mimophis mahfalensis, Oxyrhabdium leporinum, Pareas carinatus, Psammodynastes pulverulentus, Pseudoxenodon macrops, Pseudoxyrhopus heterurus, Sibynophis collaris, Stegonotus admiraltiensis, Toxicocalamus goodenoughensis, Trimeresurus albolabris, and Tropidonophis doriae. From these new genome assemblies, we extracted thousands of loci commonly used in systematic and phylogenomic studies on snakes, including target-capture datasets composed of UCEs and AHEs, as well as traditional Sanger loci. Phylogenies inferred from the two target-capture loci datasets were identical with each other, and strongly congruent with previously published snake phylogenies. To show additional utility of these non-model genomes for investigative evolutionary research, we mined the genome assemblies of two New Guinea island endemics in our dataset (Stegonotus admiraltiensis and Tropidonophis doriae) for the ATP1a3 gene, a thoroughly researched indicator of resistance to toad toxin ingestion by squamates. We find that both these snakes possess the genotype for toad toxin resistance despite their endemism to New Guinea, a region absent of any toads until the human-mediated introduction of Cane Toads in the 1930s. These species possess identical substitutions that suggest the same bufotoxin resistance as their Australian congenerics (Stegonotus cucullatus and Tropidonophis mairii) which forage on invasive Cane Toads. Herein, we show the utility of short-read high coverage genomes, as well as improving the deficit of available squamate genomes with associated voucher specimens.

Genomic resources of two aquatic Lepidoptera, Elophila obliteralis and Hyposmocoma kahamanoa, reveal similarities with Trichoptera in amino acid composition of major silk genes

Abstract While most species of butterflies and moths (Lepidoptera) have entirely terrestrial life histories, ∼0.5% of the described species are known to have an aquatic larval stage. Larvae of aquatic Lepidoptera are similar to caddisflies (Trichoptera) in that they use silk to anchor themselves to underwater substrates or to build protective cases. However, the physical properties and genetic elements of silks in aquatic Lepidoptera remain unstudied, as most research on lepidopteran silk has focused on the commercially important silkworm, Bombyx mori. Here, we provide high-quality PacBio HiFi genome assemblies of two distantly-related aquatic Lepidoptera species (Elophila obliteralis (Pyraloidea: Crambidae) and Hyposmocoma kahamanoa (Gelechioidea: Cosmopterigidae)). As a step toward understanding the evolution of underwater silk in aquatic Lepidoptera, we used our two genome assemblies and compared them to published genetic data of aquatic and terrestrial Lepidoptera. Sequences of the primary silk protein, h-fibroin in aquatic moths have conserved termini and share a basic motif structure with terrestrial Lepidoptera. However, these sequences were similar to aquatic Trichoptera in that the percentage of positively and negatively charged amino acids was much higher than in terrestrial Lepidoptera, indicating a possible adaptation of silks to aquatic environments.

Cutibacterium acnes invades prostate epithelial cells to induce BRCAness as a possible pathogen of prostate cancer

AbstractBackgroundAbundant evidence suggests that chronic inflammation is linked to prostate cancer and that infection is a possible cause of prostate cancer.MethodsTo identify microbiota or pathogens associated with prostate cancer, we investigated the transcriptomes of 20 human prostate cancer tissues. We performed de novo assembly of nonhuman sequences from RNA‐seq data.ResultsWe identified four bacteria as candidate microbiota in the prostate, including Moraxella osloensis, Uncultured chroococcidiopsis, Cutibacterium acnes, and Micrococcus luteus. Among these, C. acnes was detected in 19 of 20 prostate cancer tissue samples by immunohistochemistry. We then analyzed the gene expression profiles of prostate epithelial cells infected in vitro with C. acnes and found significant changes in homologous recombination (HR) and the Fanconi anemia pathway. Notably, electron microscopy demonstrated that C. acnes invaded prostate epithelial cells and localized in perinuclear vesicles, whereas analysis of γH2AX foci and HR assays demonstrated impaired HR repair. In particular, BRCA2 was significantly downregulated in C. acnes–infected cells.ConclusionsThese findings suggest that C. acnes infection in the prostate could lead to HR deficiency (BRCAness) which promotes DNA double‐strand breaks, thereby increasing the risk of cancer development.

Whole genome sequence resources for Pseudomonas amygdali pv. sesami and Xanthomonas arboricola isolated from sesame (Sesamum indicum) in North Carolina in 2022

Sesame (Sesamum indicum) is a specialty with increasing interest as an alternative, rotational crop for US agricultural systems. This crop is susceptible to several pathogens, including Pseudomonas and Xanthomonas species. Two bacterial isolates were recently cultured from sesame in North Carolina symptomatic of bacterial spot and identified as Pseudomonas amygdali pv. sesami (Pas) and Xanthomonas sp. These isolates were used for hybrid whole genome sequencing and assembly using Illumina and PacBio to develop robust genomic resources for these disease-causing strains. The results provide the most complete Pas genome available and placed the Xanthomoas sp. isolate into the X. arboricola species designation. Improved genomic resources for pathogens of sesame are needed to accurately detect, characterize, and employ timely management of the disease.

The genome sequence of a tachinid fly, Linnaemya tessellans (Robineau-Desvoidy, 1830)

We present a genome assembly from one male Linnaemya tessellans (tachinid fly; Arthropoda; Insecta; Diptera; Tachinidae). The genome sequence is 709.9 megabases in span. Most of the assembly is scaffolded into 7 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.24 kilobases in length.

The genome sequence of the Brown Ash Ermine moth, Zelleria hepariella Stainton, 1849

We present a genome assembly from a male Zelleria hepariella (the Brown Ash Ermine; Arthropoda; Insecta; Lepidoptera; Yponomeutidae). The genome sequence is 428.8 megabases in span. Most of the assembly is scaffolded into 19 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 16.31 kilobases in length. Gene annotation of this assembly on Ensembl identified 15,718 protein coding genes.